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Abstract 

More than twenty years ago, it was predicted that nuclei can adopt interesting 
shapes, such as rods or slabs, etc., in the cores of supernovae and the crusts of neutron 
stars. These non-spherical nuclei are referred to as nuclear "pasta." 

In recent years, we have been studying the dynamics of the pasta phases using a 
method called quantum molecular dynamics (QMD) and have opened up a new as- 
pect of study for this system. Our findings include: dynamical formation of the pasta 
phases by cooling down the hot uniform nuclear matter; phase diagrams in the density 
versus temperature plane; structural transitions between the pasta phases induced by 
compression and elucidation of the mechanism by which they proceed. In the present 
article, we give an overview of the basic physics and astrophysics of the pasta phases 
and review our works for readers in other fields. 

1 Introduction 

Massive stars with masses M = 8 — (30 — 50)M Q (M & is the mass of the sun) end their 
lives in the most spectacular and violent events in the Universe: explosions called super- 
novae, which are driven by gravitational collapse of the stellar core (see, e.g., Ref. l57l ). 
The mechanism for the collapse-driven supernova explosion has been a central mystery in 
astrophysics for almost half a century lfT2"1 . 

Great efforts have been made to unravel the mystery and a qualitative picture has been 
obtained so far (see, e.g., Refs. (6l |25l |60l |7l] |72l [731). We now know that after super- 
nova explosions the collapsing core becomes a dense compact object, neutron star, for 
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M = 8 — 30 M & or further collapses to a black hole for M > 30M©. However, no calcu- 
lations incorporating reasonable physical input succeeded in reproducing an explosion; in 
the calculations, shock waves stall on their way out of the collapsing core and the core does 
not explode but contracts into a black hole. These earlier studies suggest that interactions 
between neutrinos from the central region of the star and the shocked matter are important 
for the success of an explosion. This will hopefully revive the shock wave, and lead to a 
successful explosion. 

In equilibrium dense matter in supernova cores and neutron stars, the existence of nuclei 
with rod-like and slab-like shape is predicted Il20ll53l . Such nuclei with exotic shapes are 
referred to as nuclear "pasta." The nuclear pasta itself is an interesting subject from the 
point of view of materials science of dense stellar matter and it has been studied by nuclear 
physicists in addition to astrophysicists for more than twenty years. Furthermore, the pasta 
phases have recently begun to attract the attention of many researchers (see, e.g., Refs. 
Il9l l36ll and references therein). As has been pointed out in Refs. Il56ll65ll67l and elaborated 
in Refs. l22l l23l l58l . the existence of the pasta phases modifies the interaction between 
neutrinos and matter significantly. Our recent work [ 66 ] strongly suggests the possibility 
of dynamical formation of the pasta phases in collapsing cores from a crystalline lattice of 
spherical nuclei; effects of the pasta phases on the supernova explosions should be seriously 
discussed in the near future. 

In the present article, we provide an overview of the physics and astrophysical back- 
ground of the pasta phases for researchers in other fields, especially in soft condensed mat- 
ter physics. We try to show that the pasta phases can be an interesting system for many 
researchers in various fields, not only for nuclear astrophysicists. The plan of this paper is 
as follows. In the remaining part of the present section, we first give a brief explanation of 
collapse-driven supernovae, neutron stars and materials in these objects. We then describe 
the basic physics and astrophysical consequences of the pasta phases. In Section |2j we dis- 
cuss the similarity between the pasta phases and soft condensed matter. We then explain a 
theoretical framework used in our studies in Section [3 and show the results in Section [4] In 
this article, we generally set the Boltzmann constant fee = 1. 

1.1 Collapse-Driven Supernovae 

Stars evolve by burning light elements into heavier ones (see, e.g., Ref. [72]); here burning 
means nuclear fusion. The ultimate fate of the star is basically determined by its mass 
in the main sequence period, when the star is supported by thermal pressure of burning 
hydrogen nuclei (protons). Nuclear reactions of heavier nuclei require higher temperatures 
because of their greater Coulomb barrier. In the core of massive stars with M > 8M Q , 
the temperature reaches a threshold to produce iron. Elements heavier than iron cannot 
be produced by nuclear fusion since the binding energy per nucleon of iron is the greatest 
among all the elements. Therefore iron nuclei, the major final product of a chain of the 
thermally driven nuclear reactions, accumulate in the core of massive stars at the end of 
their evolution. 

As the nuclear burning proceeds and the mass of the iron core increases, two processes 
which tend to make the core unstable to collapse occur — electron capture, and photodis- 
sociation of heavy nuclei. In the iron core the density is so high (~ 4 x 10 9 g cm~ 3 ) that the 
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core is supported mainly by electron degeneracy pressure; degenerate electrons have ex- 
tremely large momenta even at zero temperature due to the Pauli exclusion principle. When 
the electron Fermi energy exceeds 3.7 MeV, electron capture on iron nuclei occurs: 

56 Fe + e~ -> 56 Mn + v e . (1) 

Since this reaction decreases the electron density, the electron degeneracy pressure de- 
creases and can no longer support the iron core. This instability leads to contraction of the 
core. The other instability is triggered by photodissociation of iron nuclei, which occurs 
when the temperature of the core is > 5 x 10 9 K: 

7 + jjjjPe -> 13a + 4n - 124.4MeV . (2) 

This is an endothermic reaction; thus it decreases the gas pressure and accelerates the col- 
lapse. 

As the collapse of the core proceeds and the central density approaches the nuclear satu- 
ration density (normal nuclear density), p Q = 0.165 nucleons fm~ 3 ~ 3 x 10 14 gem -3 , the 
equation of state suddenly becomes hard because of a strong short-range repulsion between 
nucleons. Due to this hardening, the pressure becomes sufficiently high to halt the collapse, 
causing the inner region of the core to bounce. The outer region of the core continues to fall 
towards the center at supersonic velocities. Consequently, the bouncing inner core drives a 
shock wave into the infalling outer core. The initial energy of the shock is ~ 10 51 erg and is 
enough to blow off the stellar envelope, which results in an explosion of the star. However, 
the shock wave propagating through the outer core is weakened by several processes (see, 
e.g., Ref. [60]) which decrease the pressure and dissipate the energy behind the shock front. 
Consequently, the shock wave stalls in the outer core. Neutrinos emitted from the inner core 
will heat the matter behind the shock front (a process referred to as neutrino heating). If the 
neutrino heating is efficient enough, the stalled shock can be revived, reach the surface of 
the outer core, propagate beyond the core, and finally blow off the outer layer of the star; 
thereby producing a supernova explosion. The contracted core remains as a nascent neutron 
star. 

The difference in the gravitational energy between the iron core and a neutron star is 
given by 



be core 



AE — — ( - ^9is) ~ ^2I£ „ O (10 53 ) erg , (3) 

V ^Fe core -^NS / -%fS 

where G is the gravitational constant, M core ~ 1M is the mass of the core, Rpe 
O(10 3 ) km is the initial radius of the iron core and the neutron star radius i?NS is about 10 
km. Only one percent of this energy is injected into the gas blown off and 99 % is carried 
away by neutrinos Q. In the supernova explosion, neutrinos interact with matter and will 
inject kinetic energy into it. Therefore detailed investigations of interactions between neu- 
trinos and matter in supernova cores are necessary to understand the explosion mechanism 
(see, e.g., Refs. ® El HUES 



'The gravitational energy is emitted in all six flavors of neutrinos and anti-neutrinos almost equally. How- 
ever, muon and tauon neutrinos (and anti-neutrinos) are not important for supernova explosions because they 
interact with matter extremely weakly. Thus, hereafter, the term "neutrino" denotes an electron neutrino. 
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Figure 1: (Color) Schematic diagram of a Wigner-Seitz cell of neutron star matter (NSM) 
and supernova matter (SNM) where n n , n p , n c and n v are the number densities of neutrons, 
protons, electrons and neutrinos respectively. 



As the density in the core increases during the collapse, the mean free path of neutrinos 
l v decreases mainly due to the neutrino coherent scattering from nuclei via the weak neutral 
current. The amplitude of the coherent scattering is proportional to A, where A is the mass 
number of the nucleus. Therefore its cross section a is proportional to A 2 lfl7l . while that 
of the incoherent scattering is simply proportional to A. When the neutrino wavelength is 
much longer than the radius of the nucleus, the neutrino is coherently scattered by nucleons 
in the nucleus. In the collapsing iron core, a typical value of the wavelength \ v of a neutrino 
with energy E v is 

he „ /10MeV\ 

K = — ~2irx 20fm( — —j , (4) 

and that of the nuclear radius is 

/ a \ 1/3 

r N ~ 1.2,4 1/3 fm ~ 5 fm — . (5) 



\56. 

The neutrino diffusion time scale is estimated by the random-walk relation: 

R 2 

Tdiff ~ —r > ( 6 ) 
cl u 

where R is the radius of the collapsing core and c is the speed of light. The dynamical time 
scale of the core collapse is comparable to the free-fall time scale: 

1 

core 
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where p core is the density of the core. When Tdiff is greater than Td yn (the corresponding 
density region is above 10 n ~ 12 gcm~ 3 ), neutrinos cannot escape from the inner core [55]]; 
this phenomenon is called neutrino trapping. While neutrinos produced by electron capture 
remain in the inner core for the order of 10 msec, the lepton fraction Yi, the lepton number 
per nucleon \Yl = (n c + n u ) / (n n + n p ) where n e , n u , n n and n p are the number densities 
of electrons, neutrinos, neutrons and protons], stays in a range of about 0.3 - 0.4. Trapped 
neutrinos build up a degenerate sea, which suppresses electron capture. Because the time 
scale of weak interaction processes is much shorter than Td yn ~ Tdiff, matter in the inner 
core is in (5 equilibrium with trapped neutrinos. At densities p < po, matter consists of 
nuclei, electrons (to maintain charge neutrality) and neutrinos, and we shall refer to this 
simply as supernova matter (SNM) (8] [30]]. We show particle number density profiles in 
supernova matter in Fig.[T]-(b). Since degenerate electrons are relativistic, screening effects 
are negligible and the electrons are uniformly distributed [,38ll50ll63l . There also exist the 
uniform neutrino gas and neutron-rich nuclei: clusters of protons and neutrons. In a certain 
density region below the normal nuclear density, supernova matter could consist of nuclei 
with exotic shapes such as rod-like and slab-like nuclei rather than spherical ones, which 
are referred to as nuclear "pasta" (a further explanation will be given in Section [T3] ). The 
pasta nuclei affect neutrino opacity of the supernova matter (see Section [L4l . 

1.2 Neutron Stars 

Let us see further evolution of the bounced core. After bounce, the core settles into hydro- 
static equilibrium on its dynamical time scale and a protoneutron star is formed unless the 
core is not so heavy that it collapses to a black hole. The initial radius and temperature of 
the protoneutron star are ~ 100 km and ~ 10 MeV, respectively. Just after the bounce of 
the core, the proton fraction x = n v j (n n + n p ), and the lepton fraction of matter in the core 
are relatively high (both are around 0.3). As neutrinos escape from the protoneutron star, 
it cools down and electron capture, which is blocked by the degenerate neutrinos, proceeds 
(see, e.g., Ref. [74] for a recent review of cooling of neutron stars). Consequently matter 
gets neutron rich; the proton fraction x < 0.1 at around the normal nuclear density p$. At 
the same time, the radius of the object shrinks to ~ 10 km; a neutron star is formed. 

Neutron stars are dense and compact objects supported by neutron degeneracy pressure 
and nuclear forces. The mass and the radius of a typical neutron star are ~ 1.4M and 
~ 10 km, respectively. The properties of dense matter and its equation of state are relatively 
well understood in the density region below ~ pq. Thus the theoretical picture of neutron 
star structure is reasonably established in the lower density region (see, e.g., Ref. |]33ll50ll 
for reviews). 

In Fig. |2] we show a schematic cross section of a typical cold neutron star below 10 9 
K. The outermost part, the envelope, is a layer about several meters thick which consists 
of liquid 56 Fe (if there is accretion from its companion star, the envelope is covered by 
a layer of hydrogen and helium atoms). With increasing density (i.e., proceeding to the 
deeper region), weak interactions render nuclei neutron-rich via electron captures triggered 
by a large electron chemical potential. Matter becomes solid at a density p > 10 6 gcm~ 3 , 
which consists of nuclei forming a bcc Coulomb lattice neutralized by a roughly uniform 
electrons. Then, at a density of about ~ 4 x 10 n gcm~ 3 , nuclei become so neutron rich 
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Figure 2: (Color) Schematic picture of the cross section of a neutron star. 



that the last occupied neutron levels are no longer bound; neutrons begin to drip out of these 
nuclei. These "dripped" neutrons form a superfluid amid the neutron-rich nuclei. The crys- 
talline region of the star is referred to as crust, which is divided into inner and outer crust; 
in the outer crust there are no neutrons outside nuclei, while in the inner crust there are. The 
inner crust extends from the neutron drip point to the boundary with the core at a density 
p < po ~ 3 x 10 14 g cnr 3 , where nuclei disappear and the system becomes uniform nuclear 
matter. Matter in these regions is called neutron star matter (NSM) (5]|7l|40]]. A schematic 
picture of neutron star matter at subnuclear densities is shown in Fig. [U-(a). There exist 
dripped neutrons outside the neutron-rich nuclei. Degenerate electrons are uniformly dis- 
tributed for the same reason as in the case of supernova matter. Another difference between 
supernova matter and neutron star matter is the lack of trapped degenerate neutrinos. This 
makes neutron star matter more neutron-rich than supernova matter in j3 equilibrium. In the 
deepest region of the inner crust corresponding to subnuclear densities (i.e., p < po), pasta 
nuclei could appear (see Section 1 1 .3 1 for a further explanation). Although crusts of neu- 
tron stars are relatively thin (~ 1 km), they influence many observed phenomena |[33ll50l . 
Effects of the pasta phases on these phenomena will be discussed in Section fOl 

The physics of the density region above several times normal nuclear density is quite 
uncertain, and a variety of hadronic phases have been proposed, such as hyperonic matter, 
pion or kaon condensates, quark-hadron mixed phase and uniform quark matter, etc. (see, 
e.g., Ref. [21j and references therein for details). These are beyond the scope of the present 
article. 

1.3 What is Nuclear "Pasta"? 

In ordinary matter, atomic nuclei are roughly spherical. This may be understood in the 
liquid drop picture of the nucleus as being a result of the forces due to the surface tension 
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Figure 3: Nuclear "pasta." The darker regions show the liquid phase, in which protons and 
neutrons coexist (i.e., the nuclear matter region); the lighter ones the gas phase, which is al- 
most free of protons. Sequence (a)-(e) shows that of nuclear shape changes with increasing 
density. This figure is taken from Ref. P31l . 



of nuclear matter, which favors a spherical nucleus, being greater than those due to the 
electrical repulsion between protons, which tends to make the nucleus deform. When the 
density of matter approaches that of atomic nuclei, i.e., the normal nuclear density po, 
nuclei are closely packed and the effect of the electrostatic energy becomes comparable to 
that of the surface energy. Consequently, at subnuclear densities around p < po/2, the 
energetically favorable configuration is expected to have remarkable structures as shown in 
Fig- El the nuclear matter region (i.e., the liquid phase of mixture of protons and neutrons) 
is divided into periodically arranged parts of roughly spherical (a), rod-like (b) or slab-like 
(c) shape, embedded in the gas phase and in a roughly uniform electron gas. Besides, there 
can be phases in which nuclei are turned inside out, with cylindrical (d) or spherical (e) 
bubbles of the gas phase in the liquid phase. As mentioned in the previous section, these 
transformations are expected to occur in the deepest region of neutron star inner crusts and 
in the inner cores of collapsing stars just before the star rebounds. Since slabs and rods look 
like "lasagna" and "spaghetti", the phases with non-spherical nuclei are often referred to as 
"pasta" phases and such non-spherical nuclei as nuclear "pasta." Likewise, spherical nuclei 
and spherical bubbles are called "meatballs" and "Swiss cheese", respectively. 

More than twenty years ago, Ravenhall et at |[53l and Hashimoto et al. [20] inde- 
pendently pointed out that nuclei with such exotic shapes can be the most energetically 
stable due to the subtle competition between the nuclear surface and Coulomb energies as 
mentioned above. Let us here show this statement by a simple calculation using an incom- 
pressible liquid-drop model (see also Ref. [46]). 

We consider five phases depicted in Fig.[3l which consist of spherical nuclei, cylindrical 



2 The "Swiss cheese" phase was considered earlier in Ref. [5] than the other pasta phases. 
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nuclei, planar nuclei, cylindrical bubbles and spherical bubbles, respectively. Each phase is 
taken to be composed of a single species of nucleus or bubble at a given nucleon density p 
averaged over the space. The total electrostatic energy including the lattice energy, which 
is the contribution to the Coulomb energy coming from neighbouring nuclei, is estimated 
by the Wigner-Seitz approximation. 

In the present discussion, we assume that the nuclear matter region is incompressible 
and the nucleon density there is fixed at the normal nuclear density po = 0.165 fm -3 . We 
further assume that the proton fraction xn in the nuclear matter region is fixed at xn = 0.3 
and that the dripped neutrons can be ignored for simplicity. These assumptions correspond 
to the situation of supernova matter [see Fig. [1Kb)]. Using a liquid-drop model, the total 
energy density of matter E tot averaged over a single cell can be written as 



where it is the volume fraction of the nuclei or bubbles, wtmik i s tne bulk energy density of 
the nuclear matter region, w SUI { , wc+l. an d E e are the nuclear surface, total electrostatic, 
and electron energy densities, respectively. Under the assumption of neglecting the dripped 
neutrons, the volume fraction u of nuclei [bubbles] is given as u = p/po [u = 1 — (p/po)]. 
The electron number density n c , which is assumed to be uniform here, is n e = pox^u 
[n e = po^n(1 — u)]. As we consider incompressible nuclear matter with fixed proton 
fraction xn, in the expression of E tot [Eq. ([8])] only w SUI { and wq+l depend on the nuclear 
shape at a given p. Thus the energetically most stable nuclear shape can be determined by 
comparing w sur { + wq+l between each nuclear shape. 

The surface energy per unit volume w sur { is expressed as 



where a is the surface tension, tn is the radius (or the half width) of nuclei or bubbles and 
d is the dimensionality of the nuclear shape defined as d = 1 for slabs, d = 2 for cylinders 
or cylindrical holes and, d = 3 for spheres or spherical holes. We assume a to be constant 
and set a = 0.73 MeV fm" 2 ll52l . 

In the Wigner-Seitz approximation, one replaces the actual unit cell by a spherical one 
for d = 3, a cylindrical one for d = 2, and a planar one for d = 1, which are electrically 
neutral and have an equal volume to the actual cell. Using this approximation, the total 
electrostatic energy per unit volume wc+l yields 




uw hu]k + w sur{ + wc+l + E e (nuclei) , 

(1 - u)w hu i k + w suvi + wc+l + E e (bubbles) 



(8) 



w c +L = 2ir(ex N p r N ) 2 uf d (u) 



(10) 



with 



fd{u) 



d+2 



1 



d- 2 



2 



1 - 



du 1 - 2 ^ 



2 



) 



+ u 



(11) 



Minimizing ^ sur f + wc+l with respect to leads to a simple condition for all the five 
nuclear shapes considered here: 

«W = 2w c +l ■ (12) 
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Thus u> su rf + t^c+L is given as 

wWf + w c +l = PoCs9d(u) , (13) 

1 /3 

with (s = 3 (f^eV^) and g d {u) = u (j^d 2 f d (u)) 1/3 . 



0.215 0.355 0.645 0.785 




0.2 0.4 0.6 0.8 

p/p 



Figure 4: (Color) u> su rf + wc+l given by Eq. (fT3l) as a function of the nucleon num- 
ber density p. The symbols SP, C, S, CH, and SH stand for nuclear shapes, i.e., sphere, 
cylinder, slab, cylindrical hole, and spherical hole, respectively. The transition points lo- 
cate symmetrically with respect to the point of p = po/2 in the case of neglecting the bulk 
contribution. 

We have plotted in Fig. [4] the summation of the nuclear surface and total electrostatic 
energy densities, w suv { + wq+l, given by Eq. (PT3l) for each five nuclear shape. At densities 
below po, the system must have an inhomogeneous structure of some kind because matter 
is assumed to be incompressible here. The point at p = po, where the lines for all shapes 
and the line of w SUT { + u>c+L = cross each other, corresponds to uniform nuclear matter. 
It can be seen from Fig. [4] that the phase with spherical nuclei gives the lowest value of 
u^surf + ^c+L at lower densities below 0.215 po, and, with increasing density, the most 
stable shape of the nuclear matter region changes as 

sphere — > cylinder — ► slab — > cylindrical hole — > spherical hole — > uniform . (14) 

Even though the pasta phases were favored energetically, whether or not they are actu- 
ally present in the star depends on occurrence of the dynamical processes leading to their 
formation. Our recent dynamical simulations |[66ll67ll68ll69l have almost solved this prob- 
lem (see Section IU). 
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1.4 Astrophysical Consequences 

The presence of nuclear pasta would affect some astrophysical phenomena. It is noted that 
the pasta phases can occupy half of the total mass of the neutron star crusts [34 ] and 10 - 
20% of the mass of supernova cores in the later stage of the collapse ll58l . In this section 
we explain effects of the pasta phases on the neutrino trapping and the core dynamics of 
supernovae, pulsar glitches and cooling of neutron stars. 

1.4.1 Supernova Explosions 

Formation of the pasta phases in collapsing cores would influence neutrino transport in the 
cores. As mentioned in Section [TTT1 this is considered to be crucial for reproduction of su- 
pernova explosions by simulation studies. Due to neutrino-nucleus coherent scattering, the 
pasta phases have a large effect on the neutrino opacity of supernova matter ll65l l67l l68l . 
If one takes account of the pasta phases, they replace uniform matter close to a boundary 
between uniform and non-uniform phases; the neutrino opacity of the pasta phases is sig- 
nificantly larger than that of uniform nuclear matter. This enhancement is elaborated with 
quantum molecular dynamic (QMD) simulations in Refs. ll22l l23l l58l (see also Section [3] 
for the framework of the QMD). The cross section for coherent neutrino scattering is ap- 
proximated to be proportional to the static structure factor S nn (q) of neutrons, where q is 
the wave number of the momentum transfer. The static structure factor is obtained by the 
Fourier transform of a radial distribution function which represents two-point correlations. 
This quantity is remarkably increased near the wave number corresponding to the recipro- 
cal lattice vector; in the case of uniform matter this peak cannot be seen. We show in Fig. 
[5] the static structure factor S nn (q) calculated for nucleon distributions of T = 1, 2 and 3 
MeV at x = 0.3 and p = 0.175 po, where the phase with rod-like nuclei appears at low 
temperatures. The nucleon distributions are obtained by our QMD simulations, which will 
be explained in Section[4] A peak around q = 0.25 - 0.3 fm _1 at each temperature indicates 
enhancement of the neutrino opacity; at 1 MeV the cross section around the peak exceeds 
one hundred times that for uniform matter. With increasing temperature, the peak height 
drops because nuclei start to melt (see phase diagrams in Section H]). 

The existence of pasta nuclei also affects the hydrodynamics of supernova cores. Phase 
transition from a bcc lattice to uniform nuclear matter is considerably smoothed by the pres- 
ence of the pasta phases, which modifies the equation of state used in core collapse simula- 
tions [53]. In addition to the equation of state, the pasta phases also have an effect on the 
other hydrodynamic properties of supernova matter due to their different elastic properties 
from those of a crystalline lattice of spherical nuclei or uniform nuclear matter. Although 
the above effects are crucial for supernova explosions, it is not obvious how they affect the 
explosions in detail; this should be investigated by simulation studies of supernovae in the 
future. 

1.4.2 Pulsar Glitches 

Pulsars are neutron stars emitting periodic pulses, with a period equal to that of the rotation 
of the star. In more than 20 pulsars, sudden decreases of pulse periods, i.e., spinups of 
neutron star crusts, are observed. These phenomena are called "glitches." There is not yet 
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Figure 5: (Color) Static structure factor of neutrons S nn (q) calculated for nucleon distri- 
butions at x = 0.3, p = 0.175 po (rod phase at low temperatures), and at temperatures of 
lMeV (solid line), 2MeV (dash-dotted line) and 3MeV (dotted line). 



a consensus about the mechanism of glitches but a widely accepted theory is the vortex 
pinning model H. 

In neutron star crusts, dripped neutrons form a superfluid and neutron-rich nuclei, which 
may be regarded as a normal fluid, act as a pinning center for vortex lines in the neutron 
superfluid. In the vortex pinning model glitches represent vorticity jumps following catas- 
trophic unpinning of the vortex lines. Since migrations of vortex lines result in transport of 
angular momentum, this leads to spinups of the crusts. It is obvious that the existence of the 
nuclear pasta would have a large influence on pinning. However, the force needed to pin 
the vortices and the pinning rate have yet to be clarified completely even for a bcc lattice of 
spherical nuclei, rather than for non-spherical nuclei |[26l . 

There is another representative model of pulsar glitches: the starquake model [54]. Ac- 
cording to this model, a sudden starquake of the neutron star crust decreases the moment of 
inertia and hence increases the angular velocity. Nuclei in the crust form a solid, Coulomb 
lattice. The rigidity of the solid crust restores the force against deformation. When the stress 
in the crust reaches a critical value, the crust cracks. The existence of the pasta phases in- 
stead of the normal bcc lattice reduces the maximum elastic energy that can be stored in the 
crust because, as will be discussed in Section 12 the pasta phases have directions in which 
restoring force does not act like liquid crystals ll49l . Thus the pasta phases would reduce 
the strength of each glitch and shorten the time separation between glitches in this scenario. 

1.4.3 Cooling of Neutron Stars 

In the early epoch of the neutron star cooling, the main process of energy loss by neutrinos 
is the so-called URCA process: 

n — > p + e~ + 9 e and p + e~ — > n + v e . (15) 

These reactions are dominant in the case that matter in the neutron star is hot and non- 
degenerate. However, when the matter becomes cold and degenerate as in a neutron star, 
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the URCA process is strongly suppressed by the energy and momentum conservation laws. 

The presence of non-spherical nuclei would accelerate the cooling of neutron stars by 
opening the direct URCA reactions which are unlikely to occur for spherical nuclei ll34ll . 
This stems from the fact that, in non-spherical nuclei, protons have a continuous spectrum 
at the Fermi surface in the elongated directions. 

The following reactions are generally accepted as the mechanism of generating neutri- 
nos in cold and degenerate neutron stars: 

n + n^n + p + e~+P e and n + p + e~ — > n + n + i/ e , (16) 

which is referred to as the modified URCA process, and neutrino pair bremsstrahlung by 
electrons in the crust: 

e~ + (Z,A) ->e _ + (Z,A) + u + u . (17) 

In the modified URCA process, a bystander neutron absorbs the momentum and relaxes 
kinematical constraints. 

Pasta nuclei suppress the bremsstrahlung rate ll27l l5ll . If we assume that neutrino 
bremsstrahlung is dominated by the lowest reciprocal lattice vectors for which form factors 
do not vanish, the number of these vectors for each dimensional lattice leads to the neutrino 
emission rate of spherical, cylindrical and planar nuclei in the ratio of 6:3:1. If the nucleons 
in the core undergo a transition to a superfluid and/or superconducting state, the modified 
URCA process is suppressed by a factor of ~ e~ A ^ kBT , where A is the superfluid en- 
ergy gap. In this case bremsstrahlung is dominant and the neutrino emission is suppressed 
efficiently as a consequence of the presence of pasta phases. 

2 Pasta Phases as Soft Condensed Matter 

Nuclear systems can be regarded as complex fluids of nucleons. Complex fluids are the 
systems in which scales of the fluid and their constituents cannot be separated clearly; 
both scales are strongly connected with each other. Such characteristic can be seen in 
various soft condensed matter systems. For example, polymer systems are one of the typical 
complex fluids (description of polymer systems in the present section is mainly based on an 
instructive review [28 ]). Actually, similar spatial structures to nuclear pasta as shown in Fig. 
[3] can be seen in some polymer systems such as domain structures in solutions and melts 
of block copolymers. Thus it would be interesting to compare the hierarchical structure of 
nuclear systems with that of polymer systems. 

Polymers are made up of monomers connected by covalent bondings and their degrees 
of polymerization are quite large (> 100). Consequently, in polymer systems, the scale of 
constituents is as large as that of macroscopic phenomena such as phase separation of the 
domain structures. Thus polymer systems behave as complex fluids. In nuclear systems, 
the nature of a complex fluid is caused by strong and short range nuclear force. 

In Fig. [6l we show the hierarchical structures of nuclear systems and of polymer sys- 
tems. Nuclear pasta phases and domain structures of block copolymers can be classified 
in the semimacroscopic scale. Thus theoretical methods based on the mesoscopic or semi- 
macroscopic scales are practical for studying these systems. In the study of polymer sys- 
tems, the density functional theory and the Ginzburg-Landau theory correspond to such 
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Figure 6: Hierarchical structures of nuclear system and of polymer system. We refer to Ref. 
ll28l for polymer system. 



approaches. In nuclear physics, we have the Thomas-Fermi theory and the liquid-drop 
model, which are based on the scale of the nucleon density profile and of the nucleus, re- 
spectively. Actually, thermodynamic properties of the pasta phases and phase diagrams at 
zero-temperature have been investigated so far mainly by these approaches. 

However, the nature of complex fluids of nuclear systems sets limits on the validity of 
the above methods. They cannot be used for investigating phenomena in which microscopic 
degrees of freedom are important: formation process of the pasta phases, thermal fluctua- 
tions and melting of the pasta structures are typical examples. However, due to the short 
range nuclear interaction, characteristic ratio of the macroscopic or semimacroscopic scale 
in nuclear systems is much smaller than that of polymer systems (see Fig. ©. Such a small 
characteristic ratio enables us to study the pasta phases by microscopic approaches based 
on the nucleonic degrees of freedom. Our simulations are performed using one of these 
methods, which will be explained in the next section. 

Similarity between the pasta phases and some phases of liquid crystals should be also 
mentioned here. For the pasta phases with rod-like nuclei and slab-like nuclei, there are di- 
rections in which the system is translationally invariant. Consequently, any restoring forces 
do not act along such directions. As noted by Pethick and his coworkers |49ll50ll . this situa- 
tion is analogous with a liquid crystal rather than with a rigid solid. Elastic properties of the 
phases with rod-like nuclei and with slab-like ones are thus characterized by elastic con- 
stants used for the corresponding liquid crystal phases, i.e., columnar phases and smectics 
A, respectively. 

For the layered phases such as smectic A liquid crystals and the pasta phase with slab- 
like nuclei, the energy density E e \ due to displacement e of layers in their normal direction 
(taken to be parallel to the z-axis) can be written as lfT3l 

E "4(!) 2+ f< v ^ 2+o < e3 >' <i8) 
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where Vj_ shows the derivative on the x-y plane. The elastic constants B and K\ are 
associated with change of the layer spacing and with bending of the slabs. 

For the two-dimensional triangular lattice in the columnar phase of liquid crystals or in 
the pasta phase with rod-like nuclei, the energy density due to a two-dimensional displace- 
ment vector e = (e x , e y ) of rods running along the z-axis is given by |[T3Tl 



£cl-f(V ± .e) 2 + ^ 



< ). .,. dey \ 2 ( dex_ dey 



dx dy J \ dy dx 

(19) 

The elastic constant B is associated with uniform translational compression or dilation, C 
with transverse shear and A'3 with bending of the rods. 

Now let us discuss thermodynamic stability of the phase with slab-like nuclei. Accord- 
ing to the Landau-Peierls argument, strict long range positional order of layered phases will 
be destroyed at finite temperatures EH. Physically, the elastic constant B for compres- 
sion of the spacing between slab-like nuclei should be comparable to the Coulomb energy 
density: 

B ~ w c +l , (20) 

and the elastic constant for bending of slab-like nuclei K\ should have dependence of the 
surface energy density times squared of the internuclear spacing I: 

Ki ~ w sui{ I 2 ~ w c +L I 2 , (21) 

here we use the fact that w suv { and wq+l are comparable in the equilibrium state [Eq. (fl2l l. 
The above expressions are consistent with results of an explicit calculation in Ref. 11491 . 

For the one-dimensional layered lattice, the mean-square displacement can then be eval- 
uated within the harmonic approximation as [ 1 1 ] 

(e 2 ) ~ — ^=ln ( ^) ~ In ( ~) , (22) 

where L is the one-dimensional length scale of the extension of the slab. The Landau-Peierls 
instability for an infinite one-dimensional layered lattice can be seen by the logarithmic 
divergence of (e 2 ) for large L. 

Typical value of the interlayer spacing is I ~ 20 fm. As can be seen from Fig. 01 
wq+l in the supernova matter is ~ 10 _1 MeV fm~ 3 . Employing the Lindemann criterion 
(e 2 )/l 2 ~ 0.1 for melting of the lattice, limits of the extension L are given by Eq. d22l) as 
In (L/l) < 100 at T ~ 1 MeV and In (L/l) < 10 at T ~ 10 MeV. In the initial phase of 
the stellar collapse before the bounce of the core, the temperature of the supernova matter 
is the order of 1 MeV and the limit of L is even much larger than the radius of the core. 
Thus the Landau-Peierls instability does not matter in this stage. 

In the neutron star inner crusts, temperature is ~ 0.1 MeV. The proton fraction in nuclei 
is about one-third smaller than that for supernova matter; as can be seen from Eq. (fTOl . 
uic+L in neutron star matter is thus wc+l ~ 10~ 2 MeV fm~ 3 [other relevant quantities are 
comparable in the both cases]. Following the same way as in the above case, the limit of L 
is given by In (L/l) < 100, which is even much larger than the radius of neutron stars. The 
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phase with slab-like nuclei in neutron star crusts is thus stable in reality with respective to 
the Landau-Peierls instability. 

Let us then estimate the temperature at which the lattice of the rod-like nuclei melts. 
For the two-dimensional triangular lattice, the mean-square displacement is iPTTTl 



where A = y/2K%/(B + 2C). An explicit calculation in Ref. |[49l shows that values of 
the elastic constants are B ~ wq+l, C ~ wc+L and K% ~ 0.05 u>c+L l 2 - Using the 
Lindemann criterion and the above typical values for I and wq+l, one can estimate the 
melting temperature T m > 10 MeV for supernova matter and T m > 1 MeV for neutron 
star matter. In the above estimate, however, decrease of the surface tension due to the 
evaporation of nucleons and thermal broadening of the nuclear density profile is not taken 
into account. Thus the melting temperature for supernova matter would be smaller than 
the above value in the real situation. More elaborated calculations on the Landau-Peierls 
instability and of the melting temperature can be seen in Refs. Il43ll64ll65l . 

In closing the present section, let us briefly mention several phases which are not shown 
in Fig.[3] In block copolymer melts, various phases with complicated structures, e.g., gyroid 
phase, perforated lamellar phase and double diamond phase, etc., have been observed. As 
will be seen in Section |4j our simulations suggest that phases with multiply connected 
structures also appear in the nuclear systems (see also Ref. [50]). This point should be 
examined in the future study H 




3 Molecular Dynamics for Nucleons 

Since the seminal works by Ravenhall et al. [53 ] and Hashimoto et al. [20], properties of 
the pasta phases in equilibrium states have been investigated using various nuclear models. 
They include studies on phase diagrams at zero temperature ll34l l45l l47l l64l 1701 and 
at finite temperatures ll32l . These earlier works have confirmed that, for various nuclear 
models, the nuclear shape changes in the way of Eq. (fl4)) as predicted by Refs. II201I53T1 . 

In these earlier works, however, a liquid drop model or the Thomas-Fermi approxima- 
tion is used with an assumption on the nuclear shape (except for Ref. iTTOl ). Thus the phase 
diagram at subnuclear densities and the existence of the pasta phases should be examined 
without assuming the nuclear shape. It is also noted that at temperatures of several MeV, 
which are relevant to the collapsing cores, effects of thermal fluctuations on the nucleon 
distribution are significant. However, these thermal fluctuations cannot be described prop- 
erly by mean-field theories such as the Thomas-Fermi approximation used in the previous 
work 021 . 

In contrast to the equilibrium properties, dynamical or non-equilibrium aspects of the 
pasta phases had not been studied until recently except for some limited cases ll24l l50l . 
Thus it had been unclear even whether or not the pasta phases can be formed dynamically 
within the time scale of the cooling of neutron stars nor whether or not the formation of 

3 In the nuclear pasta phases, the long range Coulomb interaction plays an essential role. Thus a Helfrich- 
type simple local theory for the nuclear surface cannot be applied to the present problem. 




T 



(B + 2C)VXl 



(23) 
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Figure 7: Typical length scales in the simulated system. The grey regions show "pasta" 
nuclei with the typical values of the half width tn ~ 5-10 fm, and of the separation r c ~ 
10-15 fm. The typical value of the average nucleon density in the simulation box is about 
half of the normal nuclear density po = 0.165 fm -3 . 



the pasta phases and the transitions between them can be realized under non-equilibrium 
conditions in collapsing stellar cores. Especially, the latter problem is highly non-trivial 
because these processes are accompanied by drastic changes of nuclear shape. 

To solve the above problems, molecular dynamic approaches for nucleon many-body 
systems (a comprehensive overview is given by Ref. IHJ) are suitable. They treat the mo- 
tion of the nucleonic degrees of freedom and can describe thermal fluctuations and many- 
body correlations beyond the mean-field level. Molecular dynamics for nucleons were orig- 
inally developed for studying nuclear reactions, especially the dynamics of heavy ion col- 
lisions. Nowadays there are various versions including the fermionic molecular dynamics 
(FMD) lfl4l [TBI , the antisymmetrized molecular dynamics (AMD) fl4l and the quantum 
molecular dynamics (QMD) UlEl, etc., which differ in the form of the trial wave functions. 
They are now widely used in studies on nuclear structure as well as nuclear reactions. 

Among these models, QMD is the most practical for studying the pasta phases, which 
allows us to treat large systems consisting of several nuclei. The typical length scale r c of 
half of the inter-structure is r c ~ 10 fm and the density region of interest is around half 
of the normal nuclear density po = 0.165 fm~ 3 . The total nucleon number N necessary 
to reproduce n structures in the simulation box is N ~ po(2r c n) 3 (for slabs). It is thus 
desirable to use ~ 10 4 nucleons in order to reduce boundary effects. Such large systems 
are difficult to be handled by FMD and AMD, whose calculation costs increase as ~ A 4 , 
but are tractable for QMD, whose calculation costs increase as ~ N 2 like ordinary classical 
molecular dynamics calculations (we will see the reason of this difference later). 

It is also noted that we mainly focus on the nuclear structure from mesoscopic to macro- 
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scopic scales (see Section |2]), where the exchange effect would be less important 0. There- 
fore, it is expected that QMD is a reasonable approximation for studying the pasta phases 
even though it is less elaborate on the treatment of the exchange effect than FMD and AMD. 
Especially, at finite temperatures of several MeV, validity of QMD is ensured because the 
shell structure (see Refs. ifTOl l35l l48l for the shell effects in the pasta phases) caused from 
the exchange effect is washed out by thermal fluctuations above T ~ 3 MeV. 



3.1 Framework of QMD 

A general startpoint of constructing molecular dynamic approaches for nucleons is the time 
dependent variational principle: 







* 2 dt (V(t)\ih^-H\V(t)) 



(24) 



where the time integral in the right hand side is the action. If the variation 6 of the state 
or ty*(t) is unrestricted in the whole Hilbert space except for the fixed end points at 
time ti and t^, Eq. (l24l) gives the usual time dependent Schrodinger equation. Let us now 
restrict the space to be considered by taking some trial many-body state |<&{Q(i)}), which 
contains a set of dynamical variables (complex in general) Q[t) = {<7i(i), <72(£)> • • •}■ 

In QMD, we assume each single-particle state <fii of nucleon % is a Gaussian wave packet 
and take the centers of positions and momenta of the wave packets, Rj(i) and Pj(i), as a 
set of the dynamical variables Q(t) = {Rj(i), Pj(i); i = 1, • • • , N}; 



i(r) = <r| 



(2tt<7 2 ) 3 / 4 



cxp 



R,:) 2 



4a 2 



+ t r • P* 

n 



(25) 



where a is a parameter related to the extension of the wave packet in the coordinate space. 
Thus, a is set ~ O(l) fm for nucleons in nuclear matter at around the normal nuclear 
density. The total iV-nucleon wave function is constructed as a direct product of 

single-nucleon states defined by Eq. (l25l)Fl 



(26) 



Employing the time-dependent variational principle with the above trial state |$) yields 
the following equations of motion (QMD EOM): 



R; 



dU 

&H 
~3R 



(27) 



4 This can be seen by comparing the typical values of the exchange energy for these larger scales and of the 
energy difference between two successive phases with non-spherical nuclei. See, e.g., Ref. | 68] for detail. 

~ In FMD and AMD, \&{Q(t)}) is a Slater determinant of single-nucleon states instead of a direct product 
of them. Therefore, computational costs for calculating the expectation value of the Hamiltonian in FMD and 
AMD are larger than those in QMD by a factor of iV 2 . 
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where Ti.({Q(t)}) is the expectation value of the Hamiltonian H: 

H = (®\H\$) . 



(28) 



We note that the above semiclassical equations of motion, which incorporate the effect of 
the uncertainty principle, have the same form as the classical Hamilton's equations. 

The trial wave function of QMD given by Eq. (l26l > does not possess antisymmetric 
properties. Thus it is necessary to incorporate effects of the Pauli exclusion principle in 
some way to reproduce fermionic nature: one of the solutions is introducing the Pauli po- 
tential Vpauli- The Pauli potential, a two-body repulsive potential, which depends on the 
relative momentum as well as the relative coordinate distance, mimics the Pauli principle 
by reproducing the exchange repulsive force phenomenologically, i.e., by suppressing the 
phase space density of identical particles. 

There are several forms of the Pauli potential; the following Gaussian form is one ex- 
ample: 



where Cp, qo and po are parameters. The parameters qo and po determine the range of 
the phase space distance where the exchange repulsion acts. They are thus set so that their 
product equals to the volume element of the phase space: qo ■ po ~ h. 

However, we have to mention that the genuine effect of the Pauli principle caused by 
antisymmetrization is a iV-body correlation, not a two-body effect. Therefore, the Pauli 
potential which try to imitate the Pauli principle by two-body interaction cannot exactly 
describe quantum correlations and quantum statistical properties. Using this framework, 
we thus cannot obtain correct results for quantities sensitive to the quantum fluctuations in 
the many-body system, such as the specific heat at low temperatures (42). 

3.2 Model Hamiltonian 

In our studies on the pasta phases, we have used a QMD model Hamiltonian developed by 
Maruyama et al. 071 . This model Hamiltonian consists of four parts: 



where K is the kinetic energy; Vp au ii is the Pauli potential of the Gaussian form given by Eq. 
(1291 ); Vnuci is the nuclear interaction which contains the Skyrme potential (an effective zero- 
range nuclear interaction), symmetry energy, and momentum dependent potential; Vcoulomb 
is the Coulomb energy between protons. 

The parameters in the Pauli potential are determined by fitting the kinetic energy of the 
free Fermi gas at zero temperature. The above model Hamiltonian reproduces the proper- 
ties of nuclear matter around the normal nuclear density (i.e., the saturation density, binding 
energy, symmetry energy, and incompressibility), and of finite nuclei in the ground state, 
especially of heavier ones relevant to our interest 1291 [371 . Other properties of this Hamilto- 
nian, which are important for the pasta phases (e.g., the surface tension) are investigated in 
Ref. ll68ll . It is also confirmed that a QMD Hamiltonian close to this model provides a good 
description of nuclear reactions including the low energy region (several MeV per nucleon) 




(29) 



H — K + Vpauii + Vjxacl + Vcoulomb , 



(30) 



ED- 
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4 Simulations and Results 

Using the framework of QMD, we have solved the following two major questions Il66ll67l 
|68l|69l. 1) Whether or not the pasta phases are formed by cooling down hot uniform nuclear 
matter in a finite time scale smaller than that of the neutron star cooling? 2) Whether or not 
transitions between the pasta phases can occur by compression of matter during the collapse 
of a star? In the present section, we would like to review our these works. 

In our simulations, we treated the system which consists of neutrons, protons, and elec- 
trons in a cubic box with periodic boundary condition. The system is not magnetically 
polarized, i.e., it contains equal numbers of protons (and neutrons) with spin up and spin 
down. Relativistic degenerate electrons which ensure charge neutrality can be regarded as 
a uniform background because electron screening is negligibly small at relevant densities 
around po ll38l l50l 1631 . Consequently, one must take account of the long-range nature of 
the Coulomb interaction (see also a discussion in Ref. [62]). We calculate the Coulomb in- 
teraction by the Ewald summation procedure (3j [HI taking account of the Gaussian charge 
distribution of the proton wave packets (see Appendix A in Ref. [68 ]). 



4.1 Realization of the Pasta Phases from Hot Uniform Matter 

In Refs. i67ll68l . we showed that the pasta phases are produced from hot uniform nuclear 
matter and we studied phase diagrams at zero temperature. In these works, we first pre- 
pared a uniform hot nucleon gas at a temperature T ~ 20 MeV as an initial condition 
equilibrated for ~ 500 — 2000 fm/c in advance. To realize the ground state of matter, we 
then cooled it down slowly until the temperature got ~ 0.1 MeV or less for O(10 3 — 10 4 ) 
fm/c, keeping the nucleon number density constant. In the cooling process, we mainly used 
the frictional relaxation method (equivalent to the steepest descent method), which is given 
by the physically grounded QMD equations of motion [Eq. (1271 )1 plus small friction terms: 

■ _ on &H_ 

(3D 

Here the friction coefficients £r and £p are positive definite, which determine the relaxation 
time scale. 

The resulting typical nucleon distributions of cold matter at subnuclear densities are 
shown in Figs. [8] and [9] for proton fraction of matter x = 0.5 and 0.3, respectively. We 
see from these figures that the phases with rod-like and slab-like nuclei, cylindrical and 
spherical bubbles, in addition to the phase with spherical nuclei are reproduced. The above 
simulations show that the pasta phases can be formed dynamically from hot uniform matter 
within a time scale r ~ O(10 3 — 10 4 ) fm/c. 



4.2 Equilibrium Phase Diagrams 

As has been mentioned in Section [3l QMD provides a good description of thermal fluc- 
tuations and is suitable for investigating the pasta phases at finite temperatures. We show 
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Figure 8: (Color) Nucleon distributions of the pasta phases in cold matter at x = 0.5; (a) 
sphere phase, O.lpo (^box = 43.65 fm, N = 1372); (b) cylinder phase, 0.225po (-^box = 
38.07 fm, N = 2048); (c) slab phase, 0.4p (^box = 31.42 fm, N = 2048); (d) cylindrical 
hole phase, 0.5/9n (-^box = 29.17 fm, N = 2048) and (e) spherical hole phase, 0.6po 
(^box = 27.45 fm, N = 2048), where Lbox is the box size and N is the total number of 
nucleons. The whole simulation box is shown in this figure. The red particles represent 
protons and the green ones neutrons. Taken from Ref. ll68l . 



snapshots of nucleon distributions at T = 1, 2 and 3 MeV for two densities in Figs. [TOl and 
[TT1 These densities correspond to the phases with rod-like nuclei and with slab-like nuclei 
at T = 0. 

From these figures, we can see the following qualitative features (these features are 
common for x = 0.5 and 0.3): at T ~ 1.5 — 2 MeV (but snapshots for T ~ 1.5 MeV are 
not shown), the number of evaporated nucleons becomes significant; at T > 3 MeV, nuclei 
almost melt and the spatial distribution of nucleons are almost smoothed out. 

When we try to classify the nuclear structure systematically, the integral mean curvature 
and the Euler characteristic (a brief explanation is provided in Appendix [At are useful. Here 
we introduce their normalized quantities: the area-averaged mean curvature (H), and the 
Euler characteristic density x/V (V is the volume of the whole space). Using a combination 
of these two quantities calculated for nuclear surface 0, each pasta phase can be represented 



6 Nuclear surface generally corresponds to an isodensity surface for the threshold density pth — 0.5po in 
our simulations. 
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Figure 9: (Color) The same as Fig. [8j at x = 0.3; (a) sphere phase, O.lpo (^box = 
49.88 fm, N = 2048); (b) cylinder phase, 0.18p (^box = 41.01 fm, N = 2048); (c) 
slab phase, 0.35po (^box = 32.85 fm, N = 2048); (d) cylindrical hole phase, 0.5po 
(L box = 29.17 fm, N = 2048) and (e) spherical hole phase, 0.55p (^box = 28.26 fm, 
N = 2048). The red particles represent protons and the green ones neutrons. Taken from 
Ref. 1551 . 



uniquely, i.e., 



phase with < 



spherical nuclei 
cylindrical nuclei 

slab-like nuclei 
cylindrical bubbles 

spherical holes 



(H) > 0, x/V > 0, 
(H) > 0, x/V = 0, 
(H) = 0, x/V = 0, 
(H) < 0, x/V = 0, 
(H) < 0, x/V > 0. 



(32) 



We note that the value of x/V for the ideal pasta phases is zero except for the phase 
with spherical nuclei or spherical bubbles with positive x/V', negative x/V is not obtained 
for the pasta phases. 

The phase diagrams obtained for x = 0.5 and 0.3 are plotted in Figs.[l2]and[T3l respec- 
tively. Qualitative features of these phase diagrams are the same, but the phase separating- 
region surrounded by a dashed line is smaller for x = 0.3 than that for x = 0.5. As shown 
above, nuclear surface can be identified typically at T < 3 MeV (see the dotted lines) in 
the density range of interest. Thus the regions between the dotted line and the dashed line 
correspond to some non-uniform phase, which is however difficult to classify into a specific 
phase because the nuclear surface is hard to identify. 

In the region below the dotted lines at T < 3 MeV, where we can identify the nuclear 
surface, we have obtained the pasta phases with spherical nuclei [region (a)], rod-like nuclei 
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(a) T=1 MeV (b) T=2MeV (c) T=3MeV 




Figure 10: (Color) Nucleon distributions for x = 0.5, p = 0.225po at temperatures of 1, 
2 and 3 MeV. The total number of nucleons N = 16384 and the box size Lbox = 76.14 fm. 
The upper panels show top views along the axis of the cylindrical nuclei at T = 0, and the 
lower ones show side views. Protons are represented by the red particles, and neutrons by 
the green ones. Taken from Ref. ll69l . 



(a)T=1MeV (b) T=2MeV (c) T=3MeV 




Figure 11: (Color) Nucleon distributions for x = 0.5, p = OApo at temperatures of 1, 2 
and 3 MeV. ./V = 16384 and L^ ox = 62.84 fm. These figures are shown in the direction 
parallel to the plane of the slab-like nuclei at T = 0. Protons are represented by the red 
particles, and neutrons by the green ones. Taken from Ref. [69]. 
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Figure 12: (Color) Phase diagram of matter at x = 0.5 plotted in the p - T plane. The 
dashed and the dotted lines on the diagram show the phase separation line and the limit 
below which the nuclear surface can be identified, respectively. The dash-dotted lines are 
the phase boundaries between the different nuclear shapes. The symbols SP, C, S, CH, 
SH, U stand for nuclear shapes, i.e., sphere, cylinder, slab, cylindrical hole, spherical hole 
and uniform, respectively. The parentheses (A,B) denote an intermediate phase between 
A and B-phases with a multiply connected structure characterized by negative x/V. The 
regions (a)-(g) correspond to the nuclear shapes characterized by (H) and x/V as follows: 
(a) (H) > 0, x/V > 0; (b) (H) > 0, X /V = 0; (c) (H) > 0, X /V < 0; (d) (H) = 
0, x/V = 0; (e) (H) < 0, X /V < 0; (f) (H) < 0, X /V = 0; (g) (H) < 0, X /V > 0. 
Simulations have been carried out at points denoted by circles. Adapted from Ref. |69l 
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Figure 13: (Color) Phase diagram of matter at x = 0.3 plotted in the p-T plane. Adapted 
from Ref. g9). 
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Figure 14: (Color) Phase diagram of cold matter at x = 0.1. The density at which the 
area-averaged mean curvature of nuclear surface is zero is denoted by (H) = ((H) > at 
lower densities than this point and (H) < at higher densities). However, slab phase is not 
observed in our results even at this density. Simulations have been carried out at densities 
denoted by small circles. Adapted from Ref. [68]. 



[region (b)], slab-like nuclei [region (d)], cylindrical holes [region (f)] and spherical holes 
[region (g)]. In addition to these pasta phases, structures with negative x/V have been 
also obtained in the regions (c) and (e); matter consists of highly connected nuclear and 
bubble regions (i.e., sponge-like structure) with branching rod-like nuclei, perforated slabs 
and branching bubbles, etc. A detailed discussion on the phase diagrams is given in Ref. 
1691 . Effects of the finite box size on the density regions of the phases with negative x/V 
should be examined in the future. If these density regions do not decrease significantly even 
for larger systems, the phases with multiply connected structures would be likely. 

We have also investigated cold neutron-rich matter at x = 0.1 as a more realistic condi- 
tion for neutron star matter [68 ]. The resulting zero-temperature phase diagram is shown in 
Fig-El The pasta phase with rod-like nuclei is obtained around 0.2po- A striking feature is 
that a wide density region from ~ 0.25po to ~ 0.525po is occupied by phases with negative 
x/V. The structure of matter changes rather continuously from branching rod-like nuclei 
connected to each other [obtained in the lower density region of the intermediate phase 
denoted by (C,U)] to branching bubbles connected to each other [higher density region of 
the intermediate phase (C,U)]. However, in the present neutron-rich case, pasta phase with 
slab-like nuclei cannot be obtained as far as we have investigated. A further study with a 
longer relaxation time scale is necessary to determine whether or not the phase with slab- 
like nuclei is really prohibited in such neutron-rich matter in the present model. 

4.3 Structural Transitions between the Pasta Phases 

In Ref. ll66l . we approached the second question asked at the beginning of the present sec- 
tion. We perform QMD simulations of the compression of dense matter and have succeeded 
in simulating the transitions between rod-like and slab-like nuclei and between slab-like nu- 
clei and cylindrical bubbles. 

The initial conditions of the simulations are samples of the columnar phase (p = 
0.225po) an d of the laminar phase (p = 0.4po) with 16384 nucleons at x = 0.5 and 
T ~ 1 MeV. These are obtained in our previous work I691 , which are presented in Sec- 
tion 14.21 We then adiabatically compressed the above samples by increasing the density at 
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Figure 15: (Color) Snapshots of the transition process from the phase with rod-like nuclei 
to the phase with slab-like nuclei (the whole simulation box is shown). The red particles 
show protons and the green ones neutrons. After neighboring nuclei touch as shown by the 
circle in Fig.[T5T(2). the "compound nucleus" elongates along the arrow in Fig.[T5T(3). The 
box size is rescaled to be equal in this figure. Adapted from Ref. 



the average rate of ~ 1 .3 x 10 _5 po/ (fin / c) for the initial condition of the columnar phase and 
~7. 1 x 10 -6 po / (fm / c) for that of the laminar one. According to the typical value of the den- 
sity difference between each pasta phase, ~ 0-lpo (see Fig. [T2l . we increased the density 
to the value corresponding to the next pasta phase taking the order of 10 4 fm/c. This time 
scale is much longer than the typical time scale of the deformation and structural transition 
of nuclei (e.g., that of nuclear fission is ~ 1000 fm/c). Thus the simulated compression 
process is adiabatic with respect to the change of nuclear structure, so that the dynamics of 
the structural transition of nuclei observed in the simulations is physically meaningful, and 
is essentially independent of the compression rate, etc. Finally, we relaxed the compressed 
sample at p =0.405 po for the former case and at 0.490 po for the latter one. These final 
densities are those of the phase with slab-like nuclei and cylindrical bubbles, respectively, 
in the equilibrium phase diagram at T ~ 1 MeV (see Fig.[T2l. 

The resulting time evolution of the nucleon distribution is shown in Figs. [15] and [16] As 
can be seen from Fig. [15] the phase with slab-like nuclei is finally formed rFig.[T5T(6)l from 
the phase with rod-like nuclei [Fig. [I5kl*)l. The temperature in the final state is ~ 1.35 
MeV. We note that the transition is triggered by thermal fluctuation, not by spontaneous 
deformation of the rod-like nuclei: when the internuclear spacing becomes small enough 
and once some pair of neighboring rod-like nuclei touch due to thermal fluctuations, they 
fuse fFigs. [T3]-(2) and [13k 3)1. Then such connected pairs of rod-like nuclei further touch 
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Figure 16: (Color) The same as Fig. 1 for the transition from the phase with slab-like 
nuclei to the phase with cylindrical holes (the box size is not rescaled in this figure). After 
the slab-like nuclei begin to touch [see the circle in Fig. fT6T(2)l. the bridges first crosses 
them almost orthogonally as shown by the arrows in Fig.[T6"l-(3). Then the cylindrical holes 
are formed and they relax into a triangular lattice, as shown by the arrows in Fig. [T6l-(6). 
Adapted from Ref. ll66l . 



and fuse with neighboring nuclei in the same lattice plane like a chain reaction rFig.[T31-(4)l; 
the time scale of the each fusion process is of order 10 2 fm/c, which is much smaller than 
that of the density change. 

The tr ansition from the phase with slab-like nuclei to the phase with cylindrical holes is 
shown in Fig. [16] When the internuclear spacing decreases enough, neighboring slab-like 
nuclei touch due to the thermal fluctuation as in the above case. Once nuclei begin to touch 
rFig.ll6k2VI. bridges between the slabs are formed at many places on a time scale (of order 
10 2 fm/c) much shorter than that of the compression. After that the bridges cross the slabs 
almost orthogonally for a while [Fig. 0(3)]. Nucleons in the slabs continuously flow into 
the bridges, which become wider and merge together to form cylindrical holes. Afterwards, 
the connecting regions consisting of the merged bridges move gradually, and the cylindrical 
holes relax into a triangular lattice [Fig.[l6j-(6)]. The final temperature in this case is ~ 1.3 
MeV. 

Trajectories of the above processes on the plane of the integral mean curvature 
J dR HdA and the Euler characteristic x plotted in Fig. [FT] This figure shows that 
the above transitions proceed through a transient state with "sponge-like" structure, which 
gives negative x- As can be seen from Fig. [T7T-(a) [Fig. [T7T-(b)l. the value of the Euler 
characteristic starts to decrease from zero when the rod-like [slab-like] nuclei touch. It con- 
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Figure 17: (Color) Time evolution of J dR HdA and x during the simulations. The data 
points and the error bars show, respectively, the mean values and the standard deviations 
in the range of the threshold density (see Appendix |A"1) p t ^ = 0.3 - 0.5po> which includes 
typical values for the nuclear surface. The panel (a) is for the transition from cylindrical (C) 
to slab-like nuclei (S) and the panel (b) for the transition from slab-like nuclei to cylindrical 
holes (CH). Transient states are shown as [C,S] and [S,CH] for each transition. Adapted 
from Ref. [661. 



tinues to decrease until all of the rod-like [slab-like] nuclei are connected to others by small 
bridges at t ~ 9840 fm/c [~ 12000 fm/c]. Then the bridges merge to form slab-like nuclei 
[cylindrical holes] and the value of the Euler characteristic increases towards zero. Finally, 
the system relaxes into a layered lattice of the slab-like nuclei [a triangular lattice of the 
cylindrical holes]. Thus the whole transition process can be divided into the "connecting 
stage" and the "relaxation stage" before and after the moment at which the Euler character- 
istic is minimum; the former starts when the nuclei begin to touch and it takes ~ 3000 - 
4000 fm/c and the latter takes more than 8000 fm/c. 



4.4 Formation of the Pasta Phases 

In closing the present article, let us briefly show our recent results of a study on the for- 
mation process of pasta nuclei from spherical ones; i.e., a transition from the phase with 
spherical nuclei to that with rod-like nuclei. Time evolution of the nucleon distribution in 
the transition process is shown in Fig. [18] The initial condition of this simulation is a nearly 
perfect bcc unit cell with 409 nucleons (202 protons and 207 neutrons) at T ~ 1 MeV. We 
compressed the system in a similar way to that of the simulations explained in Section 1431 
The average rate of the density change in the present case is ~4.4x 10~ 6 po/(f m / c )- Since 
the two nuclei start to touch [see the circles in Fig.[T8T-(2)1. the transition process completes 
within ~ 2500 fm/c and the rod-like nucleus is formed. The final state fFig. [T8T-(4)1 is a 
triangular lattice of the rod-like nuclei. 

The present simulation has been performed using a small system; effects of the finite 
system size in this simulation should be examined. Detailed investigation of the transition 
using a larger system will be presented in a future publication [61 ]. 
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Figure 18: (Color) Snapshots of the transition process from the bcc lattice of spherical 
nuclei to the triangular lattice of rod-like nuclei (the whole simulation box is shown). The 
red particles show protons and the green ones neutrons. The box size is rescaled to be equal 
in this figure. 



5 Summary 

In this article, we have explained the basic physics and astrophysics of non-spherical nuclei, 
nuclear pasta, in compact objects. In supernova cores and neutron stars, density can be very 
high; when the density approaches that of atomic nuclei, the Coulomb force starts to matter 
and to compete against the surface tension of nuclei. Exotic nuclear shapes result from the 
competition between these two opposite effects. The presence of the pasta phases would 
affect various astrophysical phenomena. Especially, impacts on the neutrino opacity of 
collapsing stellar cores should be examined in the near future. 

We have discussed the similarity between the pasta phases and soft condensed matter 
systems. Pasta phases with rod-like nuclei and slab-like nuclei have similar elastic prop- 
erties to those of some phases of liquid crystals: columnar phase and smectic A phase, 
respectively. With the help of this concept, thermodynamic stability of the phase with slab- 
like nuclei and melting temperature of the phase with rod-like nuclei have been discussed. 

Nuclear systems can be regarded as complex fluids of nucleons and their nature of a 
complex fluid sets limits on the validity of methods commonly used in previous studies 
of pasta phases. However, a gap between the typical length scale of the semimacroscopic 
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pasta phases and that of the nucleon is not large. Due to this small gap, we can study pasta 
phases using some microscopic approach: the quantum molecular dynamics (QMD) is one 
possibility. 

We have explained the theoretical framework of molecular dynamics calculations for 
nucleons including QMD and then given an overview of our previous work. Our QMD 
simulations have shown that the pasta phases can be formed dynamically by cooling down 
hot uniform nuclear matter in a finite time scale ~ O(10 3 — 10 4 ) fm/c, which is much 
shorter than the cooling time of a neutron star. We have also shown that transitions between 
the pasta phases can occur by compression during the collapse of a star. Our latest result 
strongly suggests that the pasta phase with rod-like nuclei can be formed by compressing a 
bcc lattice of spherical nuclei; this result will be examined in detail. 

As we have shown, nuclear pasta is an interesting system for a wide spectrum of re- 
searchers not only to nuclear astrophysicists. Physics of the pasta phases draws on nuclear 
physics, condensed matter physics, and astrophysics. The pasta phases are important for 
understanding stellar collapse and neutron star formation, which are long-standing mystery 
in the Universe. 
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A The Integral Mean Curvature and the Euler Characteristic 

The integral mean curvature and the Euler characteristic (see, e.g., Ref. [39] and references 
therein for details) are powerful tools for extracting the morphological characteristics of the 
structure of nuclear matter. Suppose there is a set of regions R, where the density is higher 
than a given threshold density p t \ x . The integral mean curvature and the Euler characteristic 
for the surface of this region dR are defined as surface integrals of the mean curvature 
H = («i + K2)/2 and the Gaussian curvature G = «i«2> respectively; i.e., f dR HdA and 
X = f dR GdA/2ir, where k\ and K2 are the principal curvatures and dA is the area element 
of the surface of R. The Euler characteristic x depends only on the topology of R and is 
expressed as 

X = (number of isolated regions) — (number of tunnels) + (number of cavities) . (33) 

Their normalized quantities, the area-averaged mean curvature (H) and the Euler charac- 
teristic density are also useful. They are defined as (H) = A^ 1 J HdA, and x/^> respec- 
tively, where V is the volume of the whole space. 
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According to Eq. (l33l . sphere and spherical holes give x > 0; infinitely long rods, in- 
finitely extending slabs and infinitely long cylindrical holes give \ = 0- Multiply connected 
structures such as "sponge" give negative %. 
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